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ABSTRACT 


Preliminary  results  of  a  theoretical  study  of  precursor 
effects  in  electromagnetic  shock  tubes  are  presented.  In 
particular,  an  examination  is  made  of  the  theories  of  fast 
luminous  fronts  which  are  observed  to  precede  shock  waves 
in  such  shock  tubes.  Errors  in  a  theory  by  Paxton  and 
Fowler  are  indicated. 
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I  INTRODUCTION 


This  technical  report  contains  some  preliminary  results  of  a  theo¬ 
retical  study  of  precursor  effects  in  electromagnetic  shock  tubes.  The 
study  was  initiated  to  supplement  the  experiment  work  on  electromagnetic 
shock  tubes  that  is  being  conducted  in  the  Electromagnetic  Sciences 
Laboratory  here  at  SRI. 

Precursor  ionization  ahead  of  the  first  shock  in  electromagnetic 
shock  tubes  has  been  widely  observed.  A  much  less  widely  observed 
phenomenon  is  that  of  a  fast  luminous  front  preceding  the  first  shock. 
The  exact  relationship  between  ' these  two  phenomena  has  yet  to  be  estab¬ 
lished.  Tli  i  s  report  deals  primarily  with  results  obtained  on  the  fast 
luminous  fronts,  but  the  theory  suggested  for  the  luminous  fronts  also 
accounts  for  the  precursor  ionization  that  results  from  the  driving 
di scha rge. 

East  luminous  fronts  have  been  observed  to  precede  the  first  shock 
wave  in  electromagnetic  shock  tubes.  Josephson  and  Hales1  noted  a 
faintly  luminous  front  in  deuterium  that  traveled  at  speeds  between  30 
and  120  cm/Vsec  at  ambient  pressures  between  0.1  and  2.5  mm  llg.  A  con¬ 
ical  shock  tube  was  driven  by  a  3.2  mE  capacitor  charged  to  21  kv.  They 
speculated  on  the  possibility  that  the  front  was  due  to  deuterons  that 
were  accelerated  to  kilovolt  energies  by  instabilities  occurring  in  dis¬ 
charge.  Med  fo  rd ,  Powell,  and  Fletcher^  also  observed  a  luminous  front 
moving  ahead  of  the  first  shock  in  deuterium.  The  speed  of  the  first 
front  was  5  cm  Msec  at  an  ambient  pressure  of  4  mm  llg.  Energy  for  their 
shock  tube  was  supplied  by  a  1 0 - ytc f  capacitor  charged  to  10  kv.  They  con¬ 
cluded  that  this  front  had  the  properties  of  a  weak  R- type  ionization 
front.  In  the  terminology  of  fluid  dynamics,  such  a  front  is  classified 
as  a  weak  detonation  front. *  howler  and  llood^  observed  fast  luminous 
waves  in  hydrogen  and  argon.  They  reported  velocities  between  60  and 
400  cm  M-sec  at  pressures  between  0.  1  and  1  mm  Hg  using  voltages  from 

V.  Josephson  end  R.  W.  Hiles,  Spsce  Technology  Lsborstories  Report  STL/TR- 6 0-  0000-  1 9 31  3 ,  1  960 

(unpu  b 1 l shed ) . 
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3  to  9  kv.  A  companion  paper  by  Paxton  and  Fowler6  presents  a  theory 
for  breakdown  wave  propagation  and  relates  such  waves  with  the  observed 
fast  luminous  fronts. 

A  criticism  of  the  treatment  by  Paxton  and  Fowler  will  be  made  in 
this  note,  and  it  will  be  hypothesized  that  breakdown  waves  can  be 
treated  as  ionization  fronts  to  first  order,  with  electrical  current 
constituting  a  second-order  effect. 


G.  W.  Paxton  and  R.  G.  Fowler,  Phya.  Rev.  128,  993  (1962  ). 
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II  CRITICISM  OF  PAXTON -FOWLER  THEORY 


Paxton  and  Fowler6  consider  the  propagation  of  luminosity  fronts 
associated  with  the  electrical  breakdown  of  a  gas.  A  breakdown  wave 
front  was  treated  as  an  electron  shock  wave.  They  presumed  that  near 
the  electrode  where  the  potential  gradient  in  the  gas  is  greatest, 
ionization  of  a  small  quantity  of  gas  occurs  and  that  the  electrons 
produced  are  given  kinetic  energy  by  the  electric  field.  The  resulting 
localized  h i gh- temperatu re  electron  gas  is  considered  to  expand  rapidly, 
thus  producing  an  electron  shock  wave  which  propagates  into  the  ambient 
gas,  partially  ionizing  the  overrun  neutral  gas  molecules.  An  incon¬ 
sistency  in  their  paper  results  from  the  statement  that  “the  energy 
necessary  for  driving  the  shock  wave  is  given  di rectlv  to  the  electrons 
in  the  shock  zone  by  the  external  field;”  actually  their  treatment 
indicates — as  they  later  state  explicitly  —  that  there  is  only  a  secondary 
dependence  of  the  propagation  speed  on  the  diiect  effect  of  the  electric 
field.  The  primary  driving  mechanism  in  their  treatment  is  the  partial 
pressure  of  h i gh - t empe ra t u re  electrons  behind  the  front. 

They  use  a  one- dimens i ona 1 ,  s teady - s t a t e,  three-fluid,  hydro- 
dynamical  model  assuming  that  the  electron  pressure  is  much  greater  than 
the  partial  pressures  of  the  other  species,  that  there  is  no  electrical 
current,  and  that  there  is  negligible  heat  flow.  The  principal  criticism 
of  their  paper  lies  in  their  treatment  of  the  zero  electrical  current 
assumption.  As  a  first-order  condition  their  assumption  seems  quite 
reasonable,  but  they  express  this  condition  as  follows: 

nv  ~  N iVi  =  0  (1 ) 

where 

n  =  Number  density  of  electrons 

v  =  Flow  velocity  of  electrons  in  the  frame  in  which 
the  front  is  at  rest 

iV.  =  Number  density  of  ions 

Yj  =  Flow  velocity  of  ions  in  the  frame  in  which  the 
front  is  at  rest. 
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Actually  lliis  equation  is  the  steady-state  expression  for  charge  conser¬ 
vation  across  the  front  (i.e.,  ahead  of  the  front  there  are  no  electrons 
or  positive  ions  and  at  the  front  the  same  number  of  free  electrons  and 
singly  ionized  positive  ions  are  generated;  thus  nv  -  N ^ t  both  ahead  of 
and  behind  the  front).  The  condition  for  zero  electrical  current  should 
be  written  in  terms  of  the  corresponding  velocities  in  the  laboratory 
frame.  Thus 


n(v  -  l/#)  -  N .  (V.  -  VQ)  -  0  (2) 

where  l0  is  the  speed  of  the  front.  Except  for  the  trivial  case  in  which 
1Q  =  0,  these  two  equations  imply  that  /V_  =  n  and  Vt  =  v.  Hence  the  ions 
and  electrons  would  travel  together,  and  a  three-fluid  model  is  not  neces¬ 
sary  for  a  zero-current  model. 

As  a  result  of  the  misinterpretation  of  Eq.  (1)  the  solution  given 
by  Paxton  and  howler  does  have  current  flowing  in  the  laboratory  frame. 

And  since  the  electron  flow  is  in  the  same  direction  as  the  front  velocity 
while  the  ions  are  virtually  stationary,  this  gives  the  unlikely  result 
that  current  flows  against  the  electric  field  for  the  case  of  a  breakdown 
wave  originating  at  a  positive  electrode.  On  the  other  hand,  a  model  of 
this  kind  may  have  merit  for  the  case  of  a  precursor  wave  traveling  into 
an  essentially  field-free  region. 


4 


Ill  FIRST-ORDER  IONIZATION -FRONT  THEORY 


An  alternative  method  that  appears  promising  for  the  treatment  of 
both  breakdown  waves  and  precursor  waves  is  to  treat  these  phenomena  as 
ionization  fronts,  as  suggested  by  Medford,  Powell  and  Fletcher.2  To 
first  order,  the  luminous  front  preceding  a  shock  wave  in  an  electro¬ 
magnetic  shock  tube  will  be  considered  to  be  the  same  as  a  breakdown 
wave  in  a  gaseous  electrical  discharge.  It  is  assumed  that  the  break¬ 
down  initiates  at  the  electrode  with  the  greatest  poteiitial  gradient  and 
that  a  localized  region  of  hot  ionized  gas  is  formed.  Ionizing  radiation 
from  this  hot  gas  is  assumed  to  be  the  primary  driving  mechanism  for  an 
ionization  front  that  moves  out  from  the  hot  gas.  It  is  further  assumed 
that  to  first  order  there  is  no  electrical  current;  hence,  a  single-fluid 
model  will  be  used. 


Consider  a  one-dimensional  picture.  Monochromatic  ionizing  radia¬ 
tion  is  coming  from  the  negative  x-direction  and  is  being  absorbed  at 
the  ionization  front.  An  ionization  front  consists  of  a  comparatively 
thin  photoabsorbing  region  situated  between  transparent  ionized  gas  be¬ 
hind  the  front  and  opaque  un-ionized  gas  ahead  of  the  front.  The  front 
moves  in  the  positive  indirection  into  the  un-ionized  gas  at  a  velocity, 
,  that  is  determined  by  the  intensity  of  the  ionizing  radiation. 

The  ratio  of  specific  heats  will  be  taken  to  be  5/l  (i.e.,  monatomic 
molecules  are  assumed),  so  the  speed  of  sound,  c0,  in  the  ambient  gas  is 
(5  pQ )  (3  P0)]/2  where  pQ  is  the  ambient  pressure  and  pQ  is  the  ambient 
density.  The  one- d i men s i ona 1  fluid-dynamic  equations  for  transfer  of 
mass,  momentum,  and  energy  in  the  rest  frame  of  the  f-ront  are 


P  i" 


1  "1 


p i  +  p 


Pf>V0 


p  +  p  F2 
"  0  c  o '  o 


(3) 

(4) 


5  Pi 

2  Pi 


5 

2  P0 


+  -  p2  +  -  Q2 

2  0  2  V 


(5) 
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where  (1  2)()2  is  l lie  excess  kinetic  energy  per  unit  muss  available  after 
ionization  of  an  atom,  and  is  defined  by 


(6) 


with  m  mass  ol  an  atom,  h  I’laiick’s  constant,  v 
ionizing  radiation,  and  ionization  energy. 


Define 

P_ \_ 

Po 


freipienry  of  the 


(7) 


1  lie n  by  Kqs.  (3),  (-!•),  and  (5)  the  equation  for  €  is 


5/'  o 

—  *  »*  *  ?2 

c2  -  5 

1? 

f  0 

L  0  J 

or,  in  terms  of  the  speed  of  sound  in  the  ambient  gas, 
,3co  '  *o  *  2  "  ,3co  +  5»'5)«  »  H'J 


0 


(8) 


(9) 


Kalin’  has  enumerated  the  possible  solutions  of  Kq .  (9)  and  related  them 

to  different  kinds  of  fronts.  The  kind  of  front  that  is  applicable  here 
is  a  weak  /1-type  ionization  front.  It  is  characterized  by  a  speed  of 
propagation  that  is  supersonic  both  with  respect  to  the  ambient  gas 
ahead  of  the  front  and  the  ionized  gas  behind  the  front.  A  necessary 
condition  thus  obtained  from  Kq.  (9)  for  the  existence  of  such  a  front  is 


»’o  >  -J  [20  *  U<?2  *  9  c 2 )  ^ ] 


(10) 


In  practice  this  condition  is  easily  met  since  ionization  of  molecules 

typically  requires  radiation  at  wavelengths  of  many  hundreds  of  Angstroms 

while  for  Q  to  be  of  comparable  magnitude  to  Vg  at  108  cm/sec,  radiation 

near  one  Angstrom  would  be  required.  Any  source  with  most  of  its  ion- 

° 

izing  radiation  at  wavelengths  greater  than  1  A  would  certainly  fulfill 
this  condition.  Actually  one  would  expect  that  Vg  >>  Q,  cg  so  that 
e  =  1  and  Vg.  Thus  in  the  laboratory  frame  the  magnitude  of  flow 
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speed  of  the  ionized  fluid  would  be  much  smaller  than  the  rate  of  advance 
of  the  ionization  front,  so  that  to  a  good  approximation  the  flow  speed 
of  the  ionized  fluid  is  zero. 

A  model  of  such  an  ionization  front,  for  which  the  shape  of  the 
ionization  density  contour  is  easily  calculable,  is  as  follows:  Consider 
a  one-dimensional  model  in  which  a  source  of  ionizing  radiation  with  a 
flux  JQ  photons/cm2  sec  is  located  at  x  =  -®.  Assume  that  only  the 
neutral  molecules  are  effective  in  the  absorption  of  radiation,  that  the 
absorption  of  each  photon  gives  rise  to  a  single  ion-electron  pair  and 
that  the  argument  of  all  dependent  variables  is  given  by  £  =  x  -  V0f  for 
t  finite.  The  equation  for  the  absorption  of  the  radiation  flux  is 
given  by 

ami 

— =  -aJ(i)nnU)  ( i  1 ) 

where  (X  is  the  absorption  coefficient  and  nJC)  is  the  number  density  of 
neutral  molecules.  Using  the  boundary  condition  J(-°°)  -  J ,,  gives 

J(l)  =  JQ  exp  [-aj^nn(2)dz  j  (12) 

Let  the  ambient  density  of  neutral  molecules  be  n„  =  n  (®)  and 

0  n 

define  £  =  0  by  the  condition  nn(0)  =  (n#/2).  Then 

ja)  +'!>„<£)  =  JQ  (13) 

and  the  speed  of  the  front  is 


(14) 


from  Eqs.  (12),  (13),  and  (14),  the  integral  equation  for  the  shape  of 
the  front  is 


n 


0 


"„(£)  +  n 0  exp 


nn (z)di 


} 


(15) 
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with  Ol  s  0,  one  obi  a  ins 


".(4) 


I  +  exp  (-ttnB£) 


for  ihe  neutral  molecules.  Thus  the  number  density  of  the  ions  and 
electrons  is  given  by 


n,(0  = 


I  +  exp  (anQ£) 


(16) 
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IV  C0INCLIIDI1NG  REMARKS 


It  appears  that  the  principal  difference  between  a  breakdown  wave 
moving  between  a  pair  of  electrodes  and  a  fast  luminous  front  moving 
away  from  the  driving  discharge  in  an  electromagnetic  shock  tube  may 
lie  in  the  presence  or  absence  of  current  flow,  respectively.  As  a 
breakdown  wave  moves  out  from  the  initiating  electrode,  the  electrode 
is  in  effect  being  extended  into  the  gas.  Tims  there  is  a  redistribu¬ 
tion  of  surface  charge  over  the  advancing  surface  of  the  ionization 
front.  This  current  flow  is  relatively  small  and  can  be  considered  as 
a  second-order  effect. 

It  has  been  assumed  in  the  foregoing  argument  that  radiation  from 
the  hot  gas  is  the  driving  mechanism.  Soft  X-ray  emission  due  to 
bombardment  of  the  initiating  electrode  can  also  provide  a  contribution. 

Medford,  Powell,  and  Fletcher2  calculate  that  an  ionizing  photon 
flux  of  3.5  y  1023  photons  cm2  sec  is  required  to  account  for  n  speed 
of  5  *  10s  cm/sec  at  a  pressure  of  4  mm  llg.  It  is  interesting  to  note 
that  this  same  flux  would  give  a  speed  of  about  l.l  *  I08  cm/sec  at 
0.1  mm  llg,  which  gives  o  rde  r -o  f -magn  i  tilde  agreement  with  the  observa¬ 
tions  of  Josephson  and  Hales,'  as  well  as  Fowler  and  Hood.5 
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